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Abstract  of  Dissertation  Presented  to  the  Graduate  Council 

in  Partial  Fulfillment  of  the  Requirements  for  the 

Degree  of  Doctor  of  Philosophy 

IMAGE  INTERPRETATION  FOR  FIELDS  PRODUCED  BY 
HIGH  FREQUENCY  LINE  CURRENTS  OVER  FINITE 
CONDUCTING  MEDIA 


By 
Bruce  Eugene  Mathews 

April  18,  1961; 

Chairman:   M.  J.  Larsen 

Major  Department:   Electrical  Engineering 

The  magnetic  field  distribution  produced  by  a 
filamentary  conductor  in  a  nonconducting  medium  carrying 
sinusoldally  time  varying  current  parallel  to  a  semi- 
infinite,  finite  conducting  medium  is  considered.   It  is 
assumed  that  the  wave  length  of  the  fields  in  the  medium 
surrounding  the  filamentary  conductor  is  large  compared 
to  the  dimensions  of  the  configuration  so  that  retardation 
can  be  neglected.   The  study  is  devoted  particularly  to 
the  high  frequency  case  where  the  field  penetration  into 
the  conducting  medium  Is  relatively  snail;  that  is, 
good  reflection  is  obtained.   If  the  medium  has  infinite 
conductivity,  reflection  Is  perfect  and  the  resulting 


field  distribution  can  be  calculated  using  a  mirror  image 
approach.   The  purpose  of  this  study  is  to  derive  a 
modified  image  approach  which  can  be  used  when  the  reflec- 
tion is  not  perfect;  that  is,  the  medium  has  finite 
conductivity.   The  analysis  is  generalized  with  respect 
to  the  magnetic  permeability  of  the  medium. 

An  equation  is  derived  for  the  vector  magnetic 
potential  produced  in  the  region  above  the  conducting 
medium.  The  resulting  expression  is  interpreted  as  that 
potential  produced  by  the  source  and  its  mirror  image 
plus  a  potential  described  by  an  integral.   An  asymptotic 
series  approximation  to  the  integral  reveals  terms  which 
can  be  interpreted  as  line  multipoles  of  increasing  order 
located  at  the  mirror  image  position.   These  multipoles 
constitute  a  modification  of  the  mirror  image  and  can  be 
used  to  describe  nonperfect  reflection. 


vi 


SECTION  I 
INTRODUCTION 

General  Problem 

A  magnetic  field  problem  which  has  many  practical 
applications  involves  an  Infinitely  long  filamentary 
conductor  in  a  nonconducting  medium  carrying  sinusoidally 
time  varying  current  parallel  to  the  surface  of  a  flat 
semi-infinite  conducting  medium.   The  configuration  is 
illustrated  in  Figure  1  where  the  current  is  assumed  to 
be  leJ*^  .   This  is  the  commonly  employed  exponential 
notation  used  to  represent  sinusoidal  time  variation  at  a 
frequency  of  u  radians  per  second,  J  being  \^T. 

?or  this  conf i?:uration  to  closely  approximate  a 
practical  problem,  it  is  not  necessary  that  the  medium 
extend  to  infinity,  but  only  to  distances  large  compared 
to  the  height  of  the  conductor  above  the  surface.   The 
solution  applies  specifically  to  a  conductor  whose  diameter 
is  very  small  compared  to  the  height.   Or,  if  a  larger 
conductor  is  involved,  the  result  can  be  considered  a 
superposition  of  many  filamentary  conductors.   In  summary, 
this  specific  problem  is  basic  to  a  large  class  of  magnetic 
field  problems. 

The  author  became  interested  in  this  problem  while 
studying  analytic  methods  of  induction  heating  work-coil 
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Figure    1.      Filanentary  Conductor   Parallel   to  Seml-Inf Inlte 
Conducting  Medium. 
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design. 1'^  Induction  heating  involves  placing  a  conductor 
configuration  (called  a  work  coll)  carrying  radio-frequency 
currents  near  the  item  to  be  heated.   The  time  varying 
magnetic  fields  induce  eddy  currents  in  the  item  and 
produce  heat.   For  the  heating  to  be  effective,  the  item 
must  be  a  fairly  good  electrical  conductor.   Practical 
applications  involve  heating  both  ferromagnetic  and  non- 
ferromacrnetic  materials. 

Designing  a  work  coil  for  a  specific  test  involves 
determining  the  conductor  configuration  which  will  give 
the  desired  heating  pattern.   The  heating  pattern  is 
determined  by  the  magnetic  field  distribution.   It  should 
be  mentioned  that  the  interest  is  in  the  induction  field 
produced  by  the  current  since  the  dimensions  of  the  coils 
are  generally  small  compared  with  the  free  space  wave 
length  of  the  frequencies  involved.   Thus,  retardation  and 
the  radiation  field  can  be  neglected  in  any  analysis. 

In  Figure  1,  the  filamentary  conductor  represents 
part  of  the  work  coil  and  the  semi-infinite  conductor 
represents  the  item  to  be  heated.   The  normal  method  of 
analyzing  this  type  configuration  is  to  first  assume  that 
the  item  is  a  perfect  conductor;  that  is,  its  electrical 
conductivity  is  infinite.   This  is  usually  considered  a 
good  approximation  since  due  to  skin  effect  the  magnetic 
field  penetrates  only  a  very  small  distance  into  the  item. 


<^   Superscripts  refer  to  numbered  bibliography. 
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This,  of  course,  depends  on  the  magnitude  of  the  frequency 
and  the  electrical  conductivity.   However,  for  most 
applications  the  penetration  is  quite  small,  since  this 
is  generally  what  is  desired  when  heating  by  induction. 
The  "perfect  conductor"  approximation  merely  assum.es  the 
penetration  is  zero  which  im.plies  that  the  surface  is  a 
perfect  reflector. 

Having  made  this  assumption,  the  field  distribution 
above  the  surface  can  be  easily  calculated  by  replacing 
the  effect  of  the  reflector  by  a  filamentary  conductor 
located  at  the  mirror  imac^e  position  of  the  work-coil 
conductor.   This  image  conductor  is  assumed  to  carry  a 
current  equal  in  maoinitude  to  the  work-coil  conductor 
but  180  de5:rees  out  of  phase.   This  configuration,  shown 
in  Figure  2,  produces  the  same  boundary  conditions  as 
the  perfect  reflector.   The  procedure  then  is  to  base 
cslculatlons  of  inductance  9nd  surface  heating  on  the 
field  distribution  produced  by  the  filamentary  conductor 
and  its  image. 

An  interesting  question  arises  as  to  the  criterion 
which  makes  this  approximation  a  valid  one.   One  would 
assumie  that  such  a  criterion  would  involve  the  relative 
magnitude  of  the  penetration  of  the  fields  into  the 
surface.   Also,  it  would  seem  reasonable  that  some  modifi- 
cation of  the  perfect  image  approach  might  be  used  to 
better  the  approximation  for  those  cases  which  do  not  meet 


Figure   2.      Image   Conductor  Produced  by   Perfect   Reflection, 
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the  criterion,   A  study  of  a  modified  image  approach  is 
the  subject  of  this  investigation. 

Previous  Studies 

A  literature  search  revealed  several  studies  in 
this  general  area.   The  reflection  of  nonperfect  conductors 
has  been  studied  by  several  inves tic^ators  interested  in 
the  wave  propagation  from  dipoles  located  above  the 
finite  conducting  earth.   Stratton  summarizes  the  early 
work  of  Sommerfeld,  Van  Der  Pol,  Weyl,  and  others  on  this 
problem.   The  analysis  showed  that  part  of  the  field 
above  the  earth  could  be  interpreted  as  that  produced  by 
a  perfect  image.   However,  there  remained  an  additional 
term  involving  an  integral  which  was  difficult  to 
interpret.  Attempts  were  made  to  approximate  the  integral 
but  no  clear  interpretation  in  terms  of  images  resulted, 

Carson^  analyzed  a  configuration  exactly  as  that 
of  Figure  1.   However,  he  limited  his  study  to  nonferro- 
magnetic  conductors  since  his  interest  was  in  the  effect 
of  the  earth  on  overhead  transmission  lines.   He  was  able 
to  obtain  approximations  to  the  field  distributions  in 
the  form  of  an  asymptotic  series  but  made  no  attempt  to 
interpret  these  results  in  terms  of  Images. 


SECTION   II 
FIELD  EQUATIONS 

Blfil'ir  Conductor  Conf ipjuration 

It    1.3    convenient   Initially  to   set   up   the    problem 
in   terp-^s   of  txvfo   lon;^    panllel   f llnrrentary  linos   with 
currents   of  equal  nan;nitude   but   180   de.':crees    out   of   phase. 
This    could   represent   a    conductor   Tnd    its    return    path 
Joined   a   Ion?   distance    avay   to   fcrm   i   closed    circuit,      A 
cross-section   representation   is    shown   in   Fi.^ure    3.      The 
lower  half    plane    is    occupied  by   a  Triaterial  with  electrical 
conductivity  op   ^"^  reafrnetic    permeability  p.^.      The    upper 
half   plane    Is   nonconductin.^  with  permeability  \i-^. 

General    Field  Relationships 

The    analysis    uses    the   concept   of  the   vector  magnetic 
potential,   T,   where    the   general   relationships   between  this 
vector  and    the   electric   and  magnetic    fields   are 

V^  A   =  tJLH,  (1) 

V   •    A   =   -txe  M,  (2) 


E  =  -w  -  M,  (3) 

^t 


-le 


jwt 


r 

b 

1_ 


M-1 


-><- 


le 


jcjt 


\l2    '    ^2, 


Figure   3.      Cross   Section  of  Blfllar   Filamentary  Line 
Currents . 
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where   H  Is    the   map:netic   field   intensity  vector,    E   is    the 
electric   field   intensity  vector,   V   is    the   scalar 
electric   potential,    and  e    is    the    dielectric    constant   of 
the  medium. 

First   note    that   due    to   the   configuration  we    have 
assumed,    only  a   z-component   of  magnetic    potential,    A    , 
will    result   since   all   of   the   current    is    in   the    z-direction. 
Also,    recalling   that  we   are   dealing  with  frequencies 
Whose   free   space  wave    length   is    very   long   compared  with 
the   dimensions    of  the   configuration,   we   can   assume   that 
the   current    in   the    lines    is    uniform  along   the   z-direction 
and   in  turn,   A^    is    not   a   function  of  z.      This   means    that 
the   divergence   of    A     is    zero   and  from  equation    (2),   iX 
is    zero.      What  we   are    really   saying    is    that    the   conductors 
are   short   enough  so   that  no   appreciable   standing  wave    is 
produced.      Now,    if  we   neglect  any  static   potential,    V 
Is    zero.      Therefore,    the   equations    (1),    (2),    and    (3) 
reduce   to 

Vx  A    =  nH,  (^) 

V   .    A   =   0,  (5) 


10 

Equations  Describing  the  Vector  Magnetic  Potential 

Taking  the  curl  of  both  sides  of  (I4.)  and  using  (5), 
we  obtain 

V^A  =  -M-V  X  H.  (7) 

Prom  Maxwell's   equations   for  a  general  homogeneous  medium, 

V  X  H  =  ol  +  e  M,  (3) 

^t 

where  the  first  term  on  the  right  represents  the  conduction 

current  density  and  the  second  term  the  displacement 

current  density. 

Now  consider  the  upper  half  plane  of  cur 

configuration.   Since  it  is  nonconducting,  the  conduction 

current  is  zero  except  at  the  line  currents.   Also,  if  we 

neglect  radiation  we  are  in  effect  neglecting  displacement 

current.   Therefore,  we  are  assuming  that  the  curl  of  H 

is  zero  except  at  the  singularities  produced  at  the  line 

currents  and  equation  (?)  for  the  nonconducting  region 

becomes 

V^J^   =  0,       (x  >  0).      (9) 

where  the  subscript  "1"  Is  used  to  denote  the  field  in 
the  upper  half  plane. 

The  medium  in  the  lower  half  plane  Is  assumed  to 
be  a  fairly  good  conductor.   This  means  that  the  conduction 
current  is  much  greater  than  the  displacement  current  at 
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the  frequency  being  considered.   Therefore,  the  second 
term  in  (8)  can  be  neglected  and  (7)  for  this  region 
becomes 

where  the  subscript  "2"  refers  to  the  lower  half  plane. 
Substituting  (6)  into  (10)  and  assuming  e-J    time  variation 
gives 

V^Ap  =  jk^Ap,        (x  <  0),   (11) 


where 


,2  ^ 


k*^  =  ut\i202'  (12) 

Source  Potential 

Consider  now  the  vector  mac^netic  potential  produced 
in  the  upper  half  plane  by  the  two  filamentary  conductors, 
neglecting  for  a  moment  the  effect  of  the  lower  half  plane. 
Equation  (9)  shows  that  the  field  distribution  is  quasi- 
steady.   That  is,  the  spatial  distribution  is  the  same  as 
that  produced  by  a  steady  current.   The  resulting  field  is 
then  merely  the  product  of  the  time  variation  of  the  current 
and  the  steady  spatial  distribution.   For  line  currents, 
the  vector  magnetic  potential  has  the  familiar  logarithmic 
distribution.   Therefore,  the  two  lines  alone  produce  a 
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z-component   source    potential,   ^q^*   S^ven  by 


ti^Ie 


jut 


'sz 


In 


i+it 


"(b-x)^   +    (y-Ha)^ 
L(b-x)2    +    (y-a)2j 


(13) 


Fourier  Expansions  of  Potential  Distributions 

The  currents  which  are  induced  in  the  lower  half 
plane  will  also  contribute  to  the  potential  In  the  upper 
half  plane.   The  resulting  potential  must  also  be  a 
solution  of  (9).   A(2;ain  notins;  that  there  is  only  a  z- 
component  of  potential  and  that  this  component  is  only  a 
function  of  x  and  y,  the  Laplaclan  of  A  In  rectangular 
coordinates  is 


:h  =  (^A^)  I^  = 


dx^ 


^y 


T 


Uk) 


where   a      is    the    unit   vector   la  the    z-dlrection.      Now    (9) 
becomes 


a^AT„         5^A. 


br.' 


^  +  l_2ll=   0. 


(x  >    0).       (15) 


^y 


A  suitable  solution  of  (15)  can  be  written  In  terms  of  a 

Fourier  expansion.   For  the  symmetry  in  our  case,  a  Fourier 
sine  expansion  is  sufficient  because  the  distribution  is 
an  odd  function  of  y.   Therefore,  the  potential  distribution 
In  the  upper  half  plane  can  be  written  as 

g3_(n)e-''"  sin(yn)  dn  +  Ag^,   (x  >  0),   (16) 
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where  the  first  term  represents  the  contribution  from  the 
current  in  the  lower  half  plane  and  the  second  term  is 
due  to  the  source  currents  and  is  given  by  Equation  (13). 
The  potential  distribution  in  the  lower  half  plane 
must  be  a  solution  of  (11),  where  the  Laplacian  of  A  is 
defined  by  (II4.).  Again  using  a  Fourier  sine  expansion, 
we  obtain  as  a  generalized  solution  of  (11) 


'2z 


=  e^"^  j   g2(n)  e^/^^^  ^  3in(yn)  dn,   (x  <  0),   (1?) 


where    the    root  with   the    positive    real   part    is    used. 

The   two   functions   gT(n)    and  g2(n)    are    determined  by 
the   boundary  conditions    at  x   =   0,      These   conditions 
require   that   the   tangential  electric   and  magnetic    field 
intensity  be    continuous    across    the   boundary.      Prom    (Ij.)    and 
(6),    this    requires    that 


1_  ^^iz   ^    1_  ^Ag^ 
\x.\     ^x         V-2     ^x 


(x   =   0). 


(18) 


^Iz  "  ^2z' 


(x  =   0). 


(19) 


Before    the   boundary  conditions    are    applied,    it    is 
necessary   to   obtain  an    integral   representation  of  equation 
(13).      From  Bateman,'+ 


00 


e"P^   3in(at)    sin(pt)   It  =   1   In 


p^-t-(a^-p)^ 
L  p2+(a-3)2j 


Re   p 


Im   (±a±3) 


(20) 
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where  Re  and  Im  indicate  "the  real  part  of"  and  the 
"Imaginary  part  of"  respectively.   Therefore,  (13)  can 
be  written 


Asz  = 


^.^leJ-^  r-(b.x)n 


3in(yn)  sin(an)  ^, 
n 


(0  <  X  <  b).     (21) 


Substituting    (21)    into    (16)    and   using    (1?)    and    (18),   we 
obtain 


^1 


M-iI      -bn 


-ngT(n)    +  — i—  e  sin(an) 


=  pl.yn2+jk2  g2(n).         (22) 


Applying   the    boundary   condition    (19)    ~,lves 

g.(n)    +  -1-  e'^"^   sin(an)    =  gp(n). 
nit  ^ 

Combining    (22)    and    (23)    and   letting  V-J^^  '  ^^  * 

2n 


,    .        ^^l^      -bn      ,     ,       X 
g. (n)    =  —=—  e  3in(an} 

nit 


_  n4|i'vrr  +  Jk^ 


-1 


(23) 


(2U) 


and    (16)    becomes 


^Iz 


2n 


-1 


_n-Hi<'vh^+Jk^ 


.-(b+x)n  s3_n(yn) 


sin(an)   l£  +  A^,. 
n  ^^ 


{2S) 
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Now  note  that  by  (20),  the  second  term  in  the 
Integral  of  (25)  is 


-U 


]^f^r^-(b-^x)n  3in(yn)  sin(an)  ^ 


=  :!ii£!i!!in 


2      2 

(bH-x)  ^^(y+a) 

L (b+x)^+(y-a)^  J 


(26) 


Comparing  v/ith    (13) »    equation    (26)    describes    the   field 
produced  by   two   filamentary  currents    located   at  x  =  -b 
and  y  =   +  a;    that    is,    the  mirror   inage   positions   of  the 
two   source   currents.      Also,    the   negative   sign   indicates 
that   the   currents    in   the    imatres    are    l80  degrees    out   of 
phase  with  the    currents    in   the   sources.      These   are 
exactly   the  mirror   imaees    produced  with  perfect    reflection. 
Denoting   this    imaae    potential   described  by    (26)    as   Aj^^, 
we   can  rewrite    (25)    as 

2niIeJ^Y"  n  -(b+x)n      ,     .      . 

^Iz   "  — ^ 1 ii-------.  e  sin(yiO 

n  °  n+u'vh^+jlc^ 

•    sin(an)    ^  +  A^^    +  A^^,  (x  >    0),  (2?) 


The  first  term  in  (2?)  describes  that  potential  which  is 
not  accounted  for  by  the  mirror  Images,   The  purpose  of 
the  next  discussion  is  to  show  that  this  remaining 
potential  can  be  interpreted  as  being  produced  by  line 
multlpoles  located  at  the  image  positions. 


SECTION  III 
DERIVATION  OF  MULTIPOLE  IMAGES 

Asymptotic  Approximation  of  Inte?,ral 

For  purposes  of  discussion,  denote  the  first  term 

In  (27)  as  A^^,. 

2\i..lei^^  r^"  -(b+x)n 

A^2  =  -i [  ^  '^^     e  ^"^^^^  sin(yn) 

n    -^0  n+p.'v^n^7jk^ 


sin(an)  ^  (28) 


The  main  task  now  is  to  find  an  interpretation  of  this 
Integral  which  will  give  insight  to  its  significance  and 
also  make  more  practical  the  calculation  of  the  field 
distributions  which  it  describes. 

After  much  investigation  involving  various  methods, 
the  most  useful  approach  appears  to  be  one  involving  an 
asymptotic  approximation  to  the  Integral,   To  form  this 
approximation,  part  of  the  integral  Is  expanded  in  a 
Taylor's  series  about  n  =  0;  namely, 


16 


17 


f(n)    = 


n+p,'vA2+ik^ 


=  f(0)    ^  f'(0)n  -.  ^"(Q^"      .   f"'(0)nf   ^ 

21  31 


=  «-J«A  -Ji_  -«-J«/2 


M,  k 


_M-'k_ 


.   e-J3«A 


•2n 

1   -  i^ 
2    J 


n 
L^x'k, 


(29) 


Note    that   f  (n)    has   branch   points    at    +^Ae^^^^^   and   poles    at 
+kiJL'e'^'^/V(l   -  M.'^)"^^^,   which  means    that    (29)    has    a   finite 
radius   of  convergence   for   finite  k.      The    intec(;ral,    on   the 
other   hand,    extend    to    infinity.      The    result   of   replacing 
f(n)    by   the   series    in    (29)    and   inte^^rating    (28)    term  by 
tenn   is    an   asymptotic   series   which  will   .^ive   a  good 
approximation   for   large   values    of  k.      Since    this    is    the 
high  frequency    case,    the   series    drives    useful   results    for 
the    problem  being    Investigated.      A   good   discussion   of 
asymptotic   approximations    is    presented  by   Jeffreys,-^ 

The    terms    resulting   from  the   expansion   of  f(n)    and 
the   subsequent    integration   can  be    Identified  as    repeated 
derivatives   of  the  mirror   image    potential,    ^]_z»    ^^ 
defined  by  equation    (26),   by  noting   that 

Sx(p)  n       Jo 


•    sin(an)   ^, 
n 


(p  =  0,    1,   2,    3---) 


(30) 
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Therefore,  equation  (28)  can  be  expanded  as 


rz 


(n'k)^    &x3 


(3i; 


To  develop  further  the  3i~nlficance  of  A^.^,  it  Is  necessary 
to  derive  neanlngful  expressions  for  the  repeated  deriv- 
atives of  Aj^g.   This  operation  is  simplified  by  using 
the  theory  of  complex  variables,   V/ith  each  iraage  position, 
we  identify  a  complex  variable  coordinate  system  as  shown 
in  Figure  I].,    and  define 


"l  ^  ^1©-'®^  =  x^  +  jy;^, 


J®2 
U2  =  r2e  ^  =  Xp  +  Jy^. 


(32) 


Equation  (26)  expressing  the  nirror  inage  potential  can 
now  be  written 

Jut 


p.-!  le 

'iz  =  -^ ^1"  ^1  -  ^^  ^2^- 

en 


Noting  that 


(33) 


In  Uj_  =  In  r^  +  je^^, 

In  U2  =  In  r2  +  ^e^*  ^341 

and  observing  that  the  two  functions  in  (3l(.)  are  analytic 
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Figure  I4..      Complex  Coordinate   Systems    at    Image    Positions. 
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(except  at  the  origin),  we  can  write 

a^P^Clnr^^) 


=  Re 


Sx2 


(p) 


=  Re 


d(P)(lnui)' 


-  duo(p)   _ 


=  (-DP-^^Cp-Di 


cos 


P©1 


=  (.l)P-l(p.l),i2i2?2^ 
r^P 


(p  =  1,  2,  3,  •••)  ,         (35) 

where  Re  indicates  that  only  the  real  part  of  the 
expression  is  used.   Since  the  partial  derivatives  with 
respect  to  x.  and  x^  equal  that  with  respect  to  x, 
equation  (31)  can  be  written 

Jut 


_  t^lle' 


rz 


Q-U/k  ^  COS    9i   _    cos    9; 
\i   k 


Q-J7i/2/  C03    29i        COS    29p 
'    (n'k)2 


+   (2-    '2n    e"-^'^"^^    /  cos   3Qi        cos   39p 


(n'k)3 


^1^ 


r23 


(36) 


Potential  Due    to   Sint^le   Line   Source 

Up  to   this    point,   we   have  been  dealing  with  two   line 
current   sources.      This  was    convenient   since    it   enabled   the 
use   of  the    integral   representation  of  equation    (21).      Ive 
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can  now  investigate  a  single  source  by  letting  the 
distance  between  the  two  conductors  become  large  compared 
to  the  distance  each  is  above  the  surface.   This  in 
effect  isolates  the  two  conductors  from  each  other. 
Assuming  a»b  and  limiting  our  attention  to  the  fields 
in  the  vicinity  of  the  source  at  y  =  a,  r2  will  become 
much  larger  then  r^  and  equation  (36)  can  be  written 


_  ^l^^ 


jut 


rz 


'e-U/h  ^  «-Jt/2 

s. cos  e,  -  _e 

PL'kr, 


1  -  -2 cos  2Q-. 

(^i'krl)2       ^ 


+  (2-|i  '^)    -2 cos  39,  - 


(n'kri)3 


(37 


This,  then,  is  the  potential,  in  addition  to  the  source 
and  its  mirror  image,  which  is  produced  in  the  upper  half 
plane  by  a  single  filamentary  current.   The  total  potential 
in  the  upper  half  plane  produced  by  a  single  filament 
located  at  x  =  b,  y  =  0  is 


M-lIe 


Hut 


1« 


27r 


mli. 


r    (li'kr^) 


^H cos  e,  -  £6 cos  2e 

2 


+  2(2-^,'2)   e-J'^A 


(  tx'kr^) 


cos  30,  - 


(li'kr^) 

,  (^  >    0),  (38) 


where  the  first  term  represents  the  source  and  its  mirror 
image  and  r,  r^,  and  0^,  are  defined  in  Figure  5. 
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Figure   5.      Single    Filamentary  Line   Current. 
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Equation  (38)  is  the  result  desired  since  It 
produces  a  useful  asynptotic  approximation  for  high 
frequencies.   Note  that  convergence  is  largely  dependent 
on  the  product  (^,'kr,  ).   Only  a  few  terms  are  necessary 
to  give  a  good  approximation  if  this  product  is  much 
greater  than  unity.   As  the  product  approaches  infinity, 
the  fields  approach  a  distribution  produced  by  perfect 
reflection,   A  useful  physical  interpretation  of  this 
convergence  requirement  can  be  obtained  by  noting  that 

where  62  ^^  ^^^  "depth  of  penetration,"  a  factor  commonly 
used  to  give  a  measure  of  the  high  frequency  skin  effect. 
Rapid  convergence,  then,  requires 

UL'kr^  »  1  (I4.O) 

or, 

r,  »  1.  lii  62  .  (UD 

Equation  (i^l)  is  interesting  since  it  shows  that  perfect 
reflection  is  not  justified  entirely  by  a  penetration 
small  compared  to  the  other  dimensions.   The  ratio  of 
magnetic  permeabilities  also  enter  into  the  criterion. 
However,  this  is  entirely  reasonable  when  one  considers 
how  the  field  distribution  is  dependent  on  this  ratio  even 
with  steady  currents. 


2h 

Line  I^ultipole  Interpretation 

To  give  a  physical  interpretation  of  the  terms 
in  equation  (38),  assume  the  conducting  medium  in  the 
lower  half  plane  is  removed  and  a  pair  of  long  fila- 
mentary conductors  spaced  a  distance  2d  apart  with 
currents  180  degrees  out  of  phase  are  placed  near  the 
image  position  as  shown  in  Figure  6,   Again,  using  the 
complex  variable  notation  of  equation  (32),  the  vector 
magnetic  potential  produced  by  the  pair  is 


A   =  ^1^1^ 
2. 


jut 


Re 


In 


u^+d 
u, -d 


_  ^^IH 


I.e^'^^ 


Re 


2n 


In 


1  + 


2d 
un+d 


(1;2) 


Now,  if  we  let  d<<r,  , 


M-l^l* 


J<jt 


Re 


2r. 


2d 
L^IJ 


jjj   HilieJ'^^  /  2d  cos  ©1  \ 


(U3) 


With  the  restriction  that  d«r,,  the  configuration  of 
Figure  6  forms  a  line  dipole  and  equation  (i|3)  describes 
the  resulting  field.   Therefore,  the  second  tenn  in 
equation  (33)  represents  such  a  line  dipole  with  a  current 
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(x.y) 


-<5)- 


Flgure  6,   Line  Dlpole  Current  at  Image  Position. 
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h  = 


ikh) 


Note    that    the    dipole    current    is   not    In  phase   v;ith  the 
source    current. 

To    interpret    the    third   term  in    (38),    consider  a 
line   quadrupole    as    shown   in   Figure   7.      Using   the   same 
approach  as  with  the   dipole,    the    resulting   potential    is 


A„    =  ne 

z  2n 


in    f       ('^l^d)(ui-d) 
(ui+jd)(ui-Jd) 


tiil2e 


j(jt 


Re 


2n 


In        1   - 


2d' 


u^2^d2 


(US) 


Again   letting   d<<  r, , 


^^^1^2®*^'^^     /      2ci2   cos   29^ 
z  ~ 

2jt 


(U6) 


The  third  term  in  OS)  then  represents  this  type  quadrupole 
with  current 


Id  = 


_  le 


-Jn/2 


(dp,  k) 


'v^2 


(1+7) 


By  a  siniliar  approach,  the  third  term  in  (35) 
represents  a  line  multipole  involving  six  currents  as 
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^Ige 


j(jt 


-l2e 


jut 


(x,y) 


Figure  7.  Line  Quadrupole  at  Image  Position. 
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shown  in  ?lc!;ure  8  with  current 


I3  =  (2-n'2)  le-J^^^   .  (1,8) 

(dn'k)3 


The  multipole  interpretation  is  more  meaningful 
if  we  let 

d  -  l/(n'k) 

=  ^^  52   .  (1;9) 

\/2  M-i 

Using  this  spacing  gives  multipole  currents  of  the  same 
order  of  magnitude  as  the  source  current.   Note  that 
(I4.9)  still  satisfies  the  requirement  d«r]_  if  the  high 
frequency  approximation  is  valid  as  stated  by  equation  {kD 

Reflected  Impedance  Produced  by  Multipoles 

Another  view  of  the  significance  of  the  asymptotic 
expansion  can  be  obtained  by  viewing  the  effect  of  each 
multipole  in  terms  of  its  contribution  to  the  electrical 
impedance  of  the  line  source.   From  equations  (6)  and 
(37),  the  electric  field  produced  by  the  multipoles, 
Ej,2{r,,©.),  at  the  source  is 
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Figure  8.   Line  Multlpole  with  Six  Currents. 
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Ej.2(2b,0)  =  - 


hk^^{2h,0) 


-jungle 


,-J.A     ^-JTr/2 


+  (2.u'2)  ^l^h/^- 


_n'k2b     (|j,'k2b)2 

3.A 
(n'k2b)^ 


(50) 


The  impedance  per  unit  length,  Z^.,  due  to  this  electric 
field  can  be  defined  in  terms  of  the  voltage  drop 
produced,  as 

Z.J.  =  -  -r^         >    '.   =  R   +  JUL   .  (51) 

Ie>*  ^ 

With  this  approach,  we  obtain 


.-lidL^l 


_n'k2bv^    (n'k2b)2   v/2(n'k2b)3 


(52) 


J'jLr.  = 


Uli]^ 


(2-M.'^) 


n'k2bv/^   ^(iJi'k2b)3 


(53) 


It  should  be  emphasized  that  the  electric  field 
produced  by  the  source  itself  and  its  mirror  image  also 
contribute  to  the  total  impedance.   This  is  the  impedance 
which  would  result  if  perfect  reflection  were  obtained. 
Using  a  coupled  circuit  analogy,  the  impedance  given  by 
equation  (5l)  can  be  considered  an  additional  "reflected 
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impedance"  due  to  the  currents  which  actually  flow  in 
the  conducting  medium.   The  resistance  term,  R  ,  is 
due  to  the  power  loss  and  the  inductance,  L^,,  is  due  to 
the  energy  stored  in  the  magnetic  field  inside  the 
conducting  medium. 

The  first  tern  in  equations  (52)  and  (53)  is 
supplied  by  the  dipole  image.   It  is  interesting  to 
note  that  this  is  exactly  the  approximation  obtained  if 
one  first  calculates  the  field  distribution  at  the 
surface  of  the  conducting  medium  based  on  perfect 
reflection,  and  then  uses  this  surface  field  distribution 
to  calculate  the  losses  and  stored  energy  in  the  medium. 
This  approach  is  often  used  in  problems  of  this  type. 

The  quadrupole  supplies  the  second  term  in  (52) 
but  does  not  contribute  to  (53).   The  multipole  involving 
six  sources  produces  an  impedance  defined  by  the  last 
term  in  (52)  and  (53) . 


SECTION  IV 
SUMMARY  AND  CONCLUSIONS 

Practical  approximations  to  the  fields  produced 
by  high  frequency  currents  and  good  conducting  media 
can  often  be  obtained  by  assuming  perfect  reflection. 
This  analysis  has  attempted  to  shed  some  light  on  the 
basis  for  this  assumption  by  considering  the  specific 
case  of  filamentary  line  currents  parallel  to  a  flat 
semi-infinite,  finite  conducting  m.edium.   If  perfect 
reflection  were  obtained,  the  effect  of  the  nedium  could 
be  replaced  by  mirror  images  of  the  line  currents.   The 
analysis  shows  that  nonperfect  reflection  can  be 
approximated  by  multlpole  images  located  at  the  mirror 
image  positions.   These  multlpoles  produce  an  asymptotic 
approximation  which  converges  fairly  rapidly  when  good 
reflection  is  obtained. 

It  would  be  interesting  to  extend  this  analysis 
to  finite  conductors  which  are  not  parallel  to  the 
surface,  other  boundary  shapes,  and  cases  Involving 
retardation  effects.   One  might  expect  to  obtain 
similiar  modified  image  approximations  for  these  cases. 
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